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We study nonlinear dynamics of exciton-polaritons in an incoherently pumped semiconductor 
microcavity with embedded weak-contrast lattice and coupled to an exciton reservoir. We elucidate 
fundamental features of non-equilibrium exciton-polariton condensate trapped in one-dimensional 
periodical potential close to zero momentum (so-called ‘Zero-state’) and to the state at the bound¬ 
ary of Brillouin zone (‘7r-state’). Within the framework of the mean-held theory, we identify dif¬ 
ferent regimes of both relaxation and oscillatory dynamics of coherent exciton-polaritons governed 
by superpositions of Bloch eigenstates within the periodic lattice. In particular, we theoretically 
demonstrate stable macroscopical oscillations, akin to nonlinear Josephson oscillations, between dif¬ 
ferent spectral components of a polariton condensate in the momenta-space. We elucidate a strong 
inhuence of the dissipative effects and the feedback induced by the inhomogeneity of incoherent 
reservoir on the dynamics of the coherent polaritons. 

PACS numbers: 71.36.+c,71.35.Gg,03.75.Kk 


I. INTRODUCTION 


Strongly correlated Bosonic particles placed in the lat¬ 
tice potentials represent indispensable tool for funda¬ 
mental studies of quantum phenomena in modern con¬ 
densed matter and solid state physic^. In atomic op¬ 
tics impressive results have been obtained due to ex¬ 
perimental achievement of Bose-Einstein condensation 
(BEC) and manipulation of ultracold atoms in optical 
lattices with different dimensionality at low enough (tens 
of microKelvins) temperature^. Nowadays atomic con¬ 
densate Josephson junctions, Mott insulator - superfluid 
quantum phase transition, topological phases and quan¬ 
tum simulation of various Bose-Hubbard type Hamiltoni¬ 
ans represent a versatile platform for quantum informa¬ 
tion technol^ies due to controllable interactions between 
the particle^. 

Last decade there was a remarkable progress in this 
area, which was related to exciton-polaritons occurring 
inside a high quality semiconductor microcavity due to 
the strong light-mater couplinJ^M^. Extremely small ef¬ 
fective masses of these composite bosons allow for the ob¬ 
servation of high teimerature non-equilibrium BEC with 
exciton-polaritons^®!. Their unique optical properties, 
such as non-parabolic dispersion and strong nonlinearity, 
are derived from the hybrid light-matter nature of these 
elementary excitations. As a result, polariton nonlinear¬ 
ity overcomes conventional optical nonlinearity in both 
strength and, most important, in speed. This is the key 
to the practical applications of exciton-polaritons, push¬ 
ing the limits of traditional photonics. Particularly im¬ 
portant nonlinear effects are optical bistabilit}^^®^ and 
parametric scattering of polaritons in planar semicon¬ 
ductor microcavities^^®^. The interplay between polari- 


tonic dispersion and strong excitonic nonlinearity allows 
for the self-localization e ffects, i.e. spontaneous forma¬ 
tion of quantized vortice^^^^, the nucleation of dark 
solitonPSHMl self-confined solutions, termed cavity 
polariton soliton^^^®^. 

In this work we study the exciton-polariton nonequilib¬ 
rium condensate confined in a weak-contrast periodical 
potential embedded into a planar microresonator driven 
by homogeneous incoherent pump [see Eig. [^a)]. The 
potential patterning of a microcavity can be achieved by 
the variety of techniques, such as reaction ion etching®!, 
mirror thickness variat ioii®!, stress applicatioiJ®!, metal 
surface depositioi]®^®!, surface acoustical wave^®!, or by 
optical mean^®!. The idea to use a periodical potential 
for the diffraction management of matter waves has been 
suggested for the c onven tional Bose-Einstein condensates 
in atomic system^®!®!. Solid-state system supporting 
exciton-polaritons is a very promising platform for the 
investigation of the spontaneous formation of coherence 
in periodical lattices. Eor instance, spontaneous build-up 
of in-phase (zero-state) and anti-phase (yr-state) “super¬ 
fluid” states have been observed in the array of weak 
periodic potential barriers in a semiconductor microcav- 
it}I®!. The formation and fragmentation of the conden¬ 
sate into the array of wires have been observed in peri¬ 
odical potentials created by surface acoustic waves®*®!. 
The existence of novel spatially localized states of coher¬ 
ent polaritonic condensates in semiconductor microcavi¬ 
ties with fabricated periodic in-plane potentials has been 
predicted®!. Note that many fundamental aspects of the 
nonlinear waves dynamics in lattices hav e bee n under¬ 
stood in purely optical dissipative system^®*®!. 

Apart from ultracold atoms forming atomic BEC the 
lifetime of low branch polaritons is in the picosecond 
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regime and comparable to the thermalization time for 
mentioned above experiments with semiconductor mi- 
crostructures. Strictly speaking, non-equilibrium open- 
dissipative behavior represents an important intrinsic fea¬ 
ture of current polaritonic solid-state systems, since po¬ 
lar itons are subjected to rapid radiative decay and their 
population maintained due to the optical pumping. The 
pumping creates a reservoir of hot excitons, which can 
relax in energy to form a polariton condensate. Even 
though polariton-polariton interactions are repulsive, the 
essentially saturable nature of exciton-polariton interac¬ 
tions may lead to effectively attractive nonlinearity for 
sufficiently low pump powerJ^. 

Various aspect of Josephson effect for the open exciton- 
photon system has been studied In particular, 

synchronized and unsynchronized phases for a polari¬ 
ton s ystem trapped in double-well potential were dis- 
cussed^^^. The influence of Josephson effect occurring 
in the nonlinear regime on polarization dynamics and 
spontaneous spatial separation of exciton-polariton con¬ 
densates with opposite polarization has been predicted^. 
At present, the main discussion is focused on the possi¬ 
bility of time enhancement of coherent oscillations occur¬ 
ring in coupled exciton-photon condensate systenP^^^. 
It has been proposed that long-lived oscillations could 
be obtained essentially due to open character of exciton- 
polariton systems and interaction with pumped reservoir. 

The objective of this paper is to make a fundamental 
theoretical investigation of exciton-polaritons in ID peri¬ 
odical potential [Fig.j^a)], taking into account their non¬ 
equilibrium character and coupling with exciton reser¬ 
voir. 

In Sec. El we describe the model of the exciton- 
polariton condensate formed in semiconductor microcav¬ 
ity in a strong coupling regime in the presence of inco¬ 
herent homogeneous optical pump. Then, we develop a 
mean-field model for three coupled harmonics. This sim¬ 
plified model proved itself as a convenient tool for obtain¬ 
ing useful analytical relations for polaritonic eigenstates 
and band-structures of the weak-contrast polaritonic lat¬ 
tices in Sec. Ill In Sec. IV we give a numerical analysis 
of different regimes of the condensate dynamics. In par¬ 
ticular, we study dynamics of condensate oscillations in 
the middle (‘Zero-state’) and at the boundary (‘tt- state’) 
of the Brillouin zone, separately. In Sec. [V| we test the 
existence and stability of the obtained nonlinear oscil¬ 
lations in the frame of a more general open-dissipative 
Gross-Pitaevskii model. 


II. THEORETICAL MODEL 

Using mean field theory and assuming the spontaneous 
formation of exciton-polariton condensate, we consider 
the open-dissipative Gross-Pitaevskii (GP) model that 
describes our incoherently pumped condensate coupled 
with an excitonic reservoir. Polariton order parameter T 
is described by GP-type equation and exciton reservoir 
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FIG. 1: (Color online) (a) Sketch of the one-dimensional 
microcavity with a periodic coating and driven by an inco¬ 
herent optical pump. A semiconductor quantum well (QW) 
is sandwiched between two Bragg mirrors (BM). (b) Band 
structure of the polariton condensate in the vicinity of the 
condensation threshold Pth = 16.335 /im“^ps“^. Eigenstates 
IJ^gr, l-J^an and figym represent the ground, antisymmetric and 
symmetric states of the polaritons with k = 0 (‘Zero-sates’), 
whereas fi- are symmetric and antisymmetric states of 
the polaritons with kl — n (‘yr-sates’). Other parameters are: 
Vt) jh — 0.1 ps“^, l = 27r//3 = 8 /xm. 
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where the reservoir induces the net gain and the blue 
shift represented by the real and imaginary parts of the 
term T{n) = ^{Rn — 7c) — igr'ri^ respectively. R = 0.01 
ps“^ /xm^ defines the condensation rate, while 7c = 0.33 
ps“^ and 7c = 0.495 ps“^ represent the decay rates of po¬ 
laritons and reservoir excitons, respectively. Constants 
Qc = 6 X 10“^ meV/xm^ and gr = 2gc characterize the 
strengths of polariton-polariton and polariton-reservoir 
interactions, respectively. In Eq. © we allow incoherent 
(non-resonant) homogeneous pumping, Pq. Here the ki¬ 
netic energy of polaritons is characterized by the effective 
mass, m, which was taken as 10“^me of the free electron 
mass. We consider a quasi-one-dimensional polariton 
condensate (for instance, confined in a microwir^^ to 
be subjected to a periodic potential approximated by the 
harmonic function V{x) = Vbcos(/3x), where [3 = 27r// 
and I is the period of modulation [see Eig. I3«)l- 
In the limit of zero potential (Vb = 0) a steady-state 
homogeneous solution (HS) of the system 0,(121 coincide 
with those known from the homogeneous cavity without 
the lattice. Eor the sake of generality we also allow HSs 
with nonzero transversal momenta 7^ 0, which have the 
form of traveling waves 




where a condensate energy is given by H/xq = 
hfi {k, |'0o|) = (fi^/2m) k^ -h 5'c|'0o|^ + gr^o. The HS be¬ 
comes nontrivial provided that the external pump com¬ 
pensates for all losses overcoming the threshold value. 
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that is Pth = "fcfriR- Coherent exciton-polariton 
and incoherent reservoir densities are given by IV^ol^ = 
(Po - Pth)he and Uro ='ydR, respectiveljEl. 

In order to illustrate the generic properties of non- 
equilibri um pola r it on dynamics we consider a simplified 
approaclP^®^ which takes into account only three main 
spatial harmonics of the weak-contrast lattice and is rel¬ 
evant to the description of three lowest energy bands. 
We expand the condensate field into a set of spatial har¬ 
monics with the common transverse wave vector /3. In¬ 
dex modulation introduces a coupling between the har¬ 
monics, which lifts the degeneracy at the intersections of 
the dispersion curves. Then, if the solution is centered 
around a momentum k the condensate wave function T 


and reservoir density n can be approximately expressed 
as 


t) « [A{t) + , (4a) 

n(a;, t) « no + (4b) 


where ^4, 5, and C are condensate amplitudes of spa¬ 
tial harmonics. The modulation of the coherent exciton- 
polaritons evokes a spatial modulation of the reservoir 
described by the terms with n± in Eqs. (4b). 


Inserting the expansions (4a) and (4b) into the 
Eqs. ([^,([^ and collecting the terms at the same har¬ 
monics, we obtain the set of mean-field equations 
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-h 'i^To^ A Qc (1^1^ + 2|5|^ + 2|C7|^) A + 2gcA*PC P 0.5Vo [P + C) + Q(^nj^P P ti—C), (5a) 
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+ iTo 5 + + 2|Ap + 2|Cp) P + QcC^A^ + 0.5Vb^ + Qn_A 


+ iTo \ C Qc (|C^P + 2|74p + 2|5p) C + g^P*A^ -|- 0.5Vo^ + Qri-^A. 


(5b) 

(5c) 


where parameter Tq = ^{Rtiq — 7c) — igr^o character¬ 
izes gain [Re(To) > 0] or damping [Re(ro) < 0] depend¬ 
ing on the homogeneous component of the reservoir no- 
The modulation of the reservoir induces an additional 
coupling between spectral components with the complex 
amplitude 0 = (yR-i-^r)- These equations have to be 
solved self-consistently with the equations for the reser¬ 
voir accounting also the spatial modulation of the reser- 


induced by the sideband components (n±): 
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Thus, a dissipative truncated model [<§ and de¬ 
scribes nonlinear evolution of the three spatial com¬ 
ponents of the exciton-polariton condensate in weak- 
contrast lattices. 


III. HOMOGENEOUS STEADY STATE 
SOLUTION 

Eirst, within the simplified model @,(11), we study 
the formation of the condensate which is just beyond 
the condensation threshold for homogeneous pump field 
(^0 > Pth)- In particular, the spatial homogeneity of the 
pump suggests that the inhomogeneous components of 
the reservoir be small no — n^o ^ n+, n_. Therefore any 
steady-stay solution satisfies approximately the following 
conservation law 

Jo = 1^(01^ + m)\^ + \c{t)\^ « ^ (7) 

with the reservoir density given by its homogeneous 
steady-state value n^o = 7c/ P- 

In the vicinity of the threshold for weak-contrast lat¬ 
tices the band gap structure is not modified so strongly 
and we can neglect the influence of exciton-exciton scat¬ 
tering forming polariton nonlinearity. In this case the 
energy (frequency) of the steady-state condensate can be 
found in the linear limit of Eqs. ^ by searching the so¬ 
lution in the form A{t) = P{t) = 

and C{t) = The solution of the resulting lin¬ 

ear eigenvalue problem provides a dispersion relation, 
i.e. dependence of the condensate energy hfioik) on the 
transversal momentum of condensate /c, see Eig. 0(b). 

An important question arises: which values of the 
transverse momentum will be chosen by the system dur- 
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ing the spontaneous build-up of the coherence? Re¬ 
cent experiments on the condensation of polaritons in 
periodical lattices show that the condensation appears 
around the high symmetry points such as zero momen¬ 
tum k = 0 point and the boundary of the Brillouin zone 
k = P /2l ^3 | 39 | 40 ] latter case the condensate is char¬ 

acterized by the phase difference between two neighbour¬ 
ing potential minima being equal to tt; this state of the 
condensate being called “7r-state”. These states, namely 
in-phase (Zero-state) and anti-phase (7r-state), reflect the 
band structure of a one dimensional polariton array and 
dynamic characteristics of meta-stable exciton-polariton 
condensates, see Fig. [^b). 

Below we focus on the detailed study of the condensate 
dynamics around the most interesting symmetry points , 
i.e the point of zero momentum k = 0 (subsec. |m^ 
and t he boundary of the Brillouin zone k = (3 /2 (subsec. 


IIIB) within the approach represented above. 


This state has a perfect sin-like shape and is spatially 
shifted at the quarter of lattice period. 

Symmetric state: The third eigenstate has a symmetric 
shape in respect to the momentum sign inversion {b = c): 


Ttllsym — + \J 4 F’o + 


a 

b 

c 


G - 4Eo/Vo 

^2G{G-4Eo/Vo) \ +1 




(13) 


We note that unlike the ground-state (8j9) the central 
component of the symmetric state converges to zero (a ^ 
0) for a vanishing potential (Vq ^ 0)- 


B. TT-states of the condensate 


A. Zero-states of the condensate 

In the limit of a vanishing potential (Vq ^0) the con¬ 
densation build-up occurs around the energy minimum, 
i.e. around the Zero-state with k = 0 [Fig. Bb)]. There¬ 
fore we start with the investigation of the condensate 
dynamics around k = 0. Solving a corresponding eigen¬ 
value problem one obtains analytical expressions for three 
eigenenergies of the condensate with the corresponding 
eigenvectors. 

Ground state: The ground eigenstate has minimal en¬ 
ergy (or frequency /i) which is given by 

= grTiro + 2-E'o “ / 2 , (8) 

where we introduce characteristic energy Eq = /Sm 
representing the “kinetic energy” of condensate at the 
boundary of the Brillouin zone, where k = /3/2. The 
corresponding eigenvector is symmetric against inversion 
of the momentum sign (b = c) and has a strong zero- 
momentum-component (|a| > |6|,|c|): 

VTo (G + ^E,/V, 

v/2G(G + 4Eo/Vo) -1 

where coefficient G = (4Eo/Vo)‘^ + 2. The normaliza- 

tion of the eigenvector’s length is given by Eq. 0. 

Antisymmetric state: The next excited eigenstate of 
the condensate is antisymmetric in respect to the mo¬ 
mentum sign inversion {b = —c): 

Qr'^rO T dF/g, (^0) 




Let us examine condensate a round t he boundaries of 
the Brillouin zone for k = /3 [Fig. [^b)]. The 

corresponding eigenvalue problem provides a cubic poly¬ 
nomial for the eigenvalue. Its exact analytical solution is 
very cumbersome. Without the loss of generality for our 
qualitative analysis we consider only two eigenstates with 
minimal eigenfrequencies. In the limit of weak-contrast 
lattices one can neglect the G component since its fre¬ 
quency substantially overcomes the frequencies of com¬ 
ponents A and B. In this approximation one finds the 
two eigenstates with smallest eigenfrequencies: 

= grriro E Eo±Vo /2 , (14) 



( 15 ) 


The sign “-h” (“—”) in Eqs. (14), represents sym¬ 
metric (antisymmetric) eigenstates relative to the cen¬ 
tral momentum k = /3 /2. In the real space both states 
have a cos-like shape but, unlike the previous cases, their 
periods are twice larger than the period of the lattice. 

As we will see below the presence of the eigenstates 
([T^,(I5) results in the persistent oscillation of the con¬ 
densate which are similar to Josephson oscillation^^. 


IV. OSCILLATION AND RELAXATION 
DYNAMICS OF THE CONDENSATE 

Superfluid properties of the exciton-polariton conden¬ 
sate in the periodic potential is significantly modified, 
taking into account non-equilibrium processeJ^. Non¬ 
linear dynamics of the coherent exciton-polaritons in the 
weak-contrast lattice is imprinted mostly by the interfer¬ 
ence and nonlinear mixing between the condensate eigen¬ 
states discussed above. The study of their stability and 
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relaxation dynamics towards the lower energetic levels 
deserves a particular interest. The dynamics of the reser¬ 
voir and its saturation play an important role here. In 
the limit of the homogeneous reservoir (no ^ n+,n_) 
the system possesses an analog of the particles conserva¬ 
tion law 0, which mimics some aspects of the conser¬ 
vative systems. The spatial inhomogeneity of the reser¬ 
voir induces additional mechanisms of the selective mode 
relaxation similar to those known from laser physics as 


spatial hole burning. In the following subsections (IV A 


and IV B) we separately study the condensate dynamics 
around the k = 0 and k = (3 /2 momenta. 


A. Zero-states of the condensate 


To illustrate the relaxation dynamics in weak-contrast 
lattices we properly prepared the initial distribution 
of exciton-polaritons in the form of the antisymmetric 
eigenstate given by Eqs. (10), (11). Experimentally it 


can be realized by launching two coherent optical beams 
with a properly chosen phase difference and tilted in op¬ 
posite directions. The initial state of the reservoir was 
taken to be spatially homogeneous, i.e. no = n^o and 
n+ = n_ =0. We also took into account the initial 
small-amplitude noise, that has been added to the ini¬ 
tial condensate profile. This noise breaks the symmetry 
of the initial eigenstate triggering the relaxation dynam¬ 
ics. Eigures[^(a) and (b) clearly demonstrate dynamics 
and eventual relaxation to the ground states [given by 
Eqs. (|^, ^]. In the ground state coherent polaritons 
gather around the minima of the periodical potential. 

We note that this relaxation dynamics is induced by 
the spatial modulation of the reservoir. Indeed, in the 
frame of the simplified model and (|^], all three po- 
lariton eigenstates experience a uniform net gain [i?e(ro)] 
provided that the reservoir modulation is negligible, i.e. 
n± = 0. Apparently the reservoir modulation {n± ^ 0) 
breaks this symmetry in favour of the ground polaritonic 
state. That is why the ground state becomes eventually 
much more populated. 

This simple relaxation dynamics changes substantially 
for lattices with the stronger potential con tras ts Vb- In¬ 
deed, the ground state described by Eqs. becomes 

unstable and leads the system to the regime of persistent 
oscillations [Eig. (c),(d)]. The oscillations appear due 
to the temporal beating between the ground st ate an d 
the excited symmetric mode population [Eqs. ( 12|13[ )], 
cf. Fig.j^b). Thus the oscillation frequency uJzero can be 
approximated as the difference in their eigenfrequencies 
lisym and (igr^ at least, in the vicinity of the condensation 
threshold: 


^zero ~ ^{4Eo/n)^ + 2{Vo/h)^. 


(16) 


Our extensive numerical analysis of the model 0 , 0 
shows that the Eq. (16) predicts accurately the oscilla¬ 
tions frequency in the vicinity of the condensation thresh¬ 
old (Po ~ Pth) where the condensate amplitudes and 




FIG. 2: (Color online) Relaxation dynamics of Zero-state 
polaritons for pumping Po = 25 /im“^ps“^ and different po¬ 
tential contrasts: (a,b) for Vo /h = 0.2 ps“^ and (c,d) for 
Volk — 0.35 ps“^. The figures (a,c) show spatial profiles of 
condensate density |^|^, whereas the figures (b,d) represent 
the amplitudes of the spatial harmonics (|A|,|P|,|C|). The 
initial state of the condensate is taken in the form of an¬ 
tisymmetric state 0 for (a,b) and in the form of ground 
state 0 for (b,d), respectively, (e) Oscillation frequencies 
vs. the pumping rate Po for different values of modulation 
depths Vq. Dark diamonds depict the analytical results given 
by Eq. (16). (f) Phase boundary of pump Po versus modu¬ 
lation depth Volk characterizing domain of persistent oscilla¬ 
tion existence (shaded area). Vertical dashed line determines 
threshold value Pth = 16.335 /im“^ps“^ of the pump. The 
lattice period is the same as in Fig. 


reservoir modulations are small [Fig.j^e)]. For a stronger 
pumping the influence of the nonlinear effects and reser¬ 
voir inhomogeneities become non-negligible. As a conse¬ 
quence the oscillations frequency increases with the pump 
amplitude Pq until the critical value where the ground 
mode of the condensate state stabilizes again and the 
oscillations disappear. This threshold of the oscillations 
appearance is determined both by potential depth Vb and 
pumping rate Pq. The domain boundary for the regime 
of persistent oscillations is shown in Fig. |^f). 

The growth of another eigenmode against the conden¬ 
sate ground state can be explained again by the selective 
saturation of the gain. Due to the spatial dependency of 
the reservoir density the gain is saturated mostly in the 
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vicinity of the condensate maxima, i.e. in the minima of 
the periodical potential. It allows another mode, whose 
peak fields are localized near the maxima of periodical 
potential [i.e. the symmetric mode given by (EHHi)], 
the opportunity to grow as well. This phenomenon is 
known from laser physics as spatial hole burning. Sim¬ 
ilar beating dynamics between two spatial modes of a 
weak-contrast lattice has been predicted for the optical 
parametric oscillatoi^^. 


B. TT-states of the condensate 


Nonlinear relaxation dynamics of the polariton con¬ 
densate can hardly be understood without considering 
TT-states of the condensate in a periodical potential, i.e. 
the states at the boundary of the Brillouin zone. Two 
tilted coherent beams with an appropriate photon energy 
(14), (15) can create required spatial harmonics with the 


momenta k = +/3/2 and k = —/d/2. Our numerical 
analysis of the model (HI,®) shows that both symmetric 
(/i+) and antisymmetric (/i_) eigenstates undergo insta¬ 
bility and develop into dynamically stable persistent os¬ 
cillations [Fig. (a),(b)]. The maxima of the condensate 
density oscillate around the bottoms of the periodical 
potential. Moreover the oscillations between two nearest 
neighboring potential sides are in phase. 


Similar to the dynamics considered above (16) the os¬ 


cillations between two TT-states, namely, symmetric /i+ 
and antisymmetric /i_ ones, exist even in the vicinity 
of the condensation threshold where the oscillation fre¬ 
quency can be approximated by a simple expression 


■ 11+ - ij.- = Vo /h. 


(17) 


Thus, coherent polaritons periodically change their col¬ 
lective momenta between two values k = +/d /2 and 
k = —yd /2 which are represented in the model ([^,® by 
the spatial harmonics A and 5, respectively [Fig. I^b)]. 
This behaviour is similar to Josephson oscillations of the 
polaritons in the double well potential welP^. However, 
unlike conventional Josephson dynamics, the oscillations 
being considered here appear in the /r-space, i.e. between 
two components with different momenta and, therefore, 
can be observed in the far-held measurements. 

Alternatively the frequency of oscillations can be 
found analytically in the limit of very weak lattice con¬ 
trasts (Vo ^ 0) where the spatial modulation of the reser¬ 
voir can be neglected {n± = 0) against the main compo¬ 
nent no = 7c/i? created by the homogeneous pump. In¬ 
deed, by setting C = 0 and n± = 0 in ([^,(|^, one reduces 
the problem to the standard system of two coupled non¬ 
linear equationJ^. In this case, supposing that during 
oscillations components A and B completely exchange 
their populations, one can obtain 


2hF{^ |m)’ 


(18) 




FIG. 3: (Color online) Relaxation dynamics of TT-states po¬ 
laritons for the potential contrasts Vo /h = 0.1 ps~^ and dif¬ 
ferent pumping rates: (a,b) for Pq = 18 /im“^ps“^ and (c,d) 
for Pq = 36 /im“^ps“^. Figures (a,c) show the spatial prohles 
of the condensate density whereas the hgures (b,d) repre¬ 

sent the amplitudes of the spatial harmonics (|A|,|H|,|C|). (e) 
Oscillation frequencies vs. the pumping rate Pq for different 
values of the modulation depths Vq. Dark diamonds depict 
the analytical results given by Eq. 0- (f) Phase boundary 
of the pump Pq versus modulation depth Vo/h. The lattice 
period is the same as in Fig. 


where F {x\m) is elliptic integral of the first order with 
m = The total density of the condensate Iq = 

\Af + \Bf can be approximated by the Eq. In the 
vicinity of the condensation threshold m 1 and, thus, 
the frequency of TT-oscillations approaches 


U) 


7T 


Ed , 

h \ 16 


-1 


(19) 


We note that for Pq Pth the oscillations become linear, 
i.e. Iq 0, and Eq. (19) takes the form Eq. (17). In fact. 


the approach given by the Eqs. ( |18| ) and (|19fis valid in 
the vicinity of threshold, since, unlike the conventional 
case of Josephson oscillations, the system is governed by 
the saturation dynamics of the reservoir rather than by 
the Kerr-like nonlinearity. 

For a stronger pumping nonlinear effects change the 
character and the periods of oscillations drastically. The 



















































7 


mean frequency decreases along with the increase 
of pump P. Simultaneously, the form of oscillations be¬ 
comes less regular, and the dynamics transforms even¬ 
tually into a quasi-periodical regime. For even stronger 
pumping (Pq > 36 /im“^ps“^ for Vq/H = 0.1 ps“^) the 
solution becomes steady-state [Fig. ic),(d)]. 

In general, the character of 7r-state oscillations and 
their frequencies are sensitive to the system parameters, 
as it has been summarized in Fig. [^e). For the lattices 
with very weak modulation contrasts (Vo/h < 0.21 ps“^) 
persistent oscillations exist in a pump interval from the 
condensation threshold (Pq « Pth) to the critical value 
where the oscillation frequency approaches zero. The ex¬ 
istence domain of the persistent oscillations are presented 
in Fig. [^f). For stronger modulation depths (Vo/h > 
0.21 ps“^), however, the threshold of 7r-oscillations shifts 
to the larger values of the pumping rates and, therefore, 
the periodical energy exchange between k = +/d /2 and 
k = —p /2 components are impossible in the linear limit. 
Notably, 7r-oscillation regime discussed above is violated 
in the upper (depicted by “multi-mode oscillations” ) 
region in Fig. if)- Being in this regime the exciton- 
polariton system exhibits oscillatory behavior for which 
the population of mode \Cf‘ becomes comparable with 
1^41^ and \B\^. 


V. NONLINEAR DYNAMICS WITHIN THE 
GROSS-PITAEVSKH DISSIPATIVE MODEL 


In the previous sections we discussed nonlinear dynam¬ 
ics in the frame of a simplified model ( 5a]|6d ). However, 
to obtain a more general picture of the relaxation and os¬ 
cillations between all available spatial components, one 
has to use the original model represented by Eqs. 0 ,( 111 ). 
In this section we prove the existence and stability of both 
types of oscillatory solutions obtained above by means of 
the extensive numerical analysis within the frame of the 
open-dissipative Gross-Pitaevskii model with the period¬ 
ical potential. 

First, we applied the appropriately tilted optical beam 
and created the initial condensate with the central mo¬ 
mentum k ^ P /2. As a result we observed persistent 
oscillations over a long time exceeding, being at least 
several tenths of nanoseconds [Fig. Qa)]. The evolution 
of the phase profile shows that coherent polaritons really 
change periodically their collective momenta between two 
values k = P/S /2 and k = —/3 /2 [Fig. (b)]. This is in 
perfect agreement with our analytical predictions given 
above [see Eq. and Eig. I^a) and (e)]. To check the 
stability of these oscillations we added some additional 
noise to the initial condensate distribution. Note that 
these oscillations are very robust and are insensible to 
the particular value of the amplitude of the initial tilted 
beam. We found that both quantitative and qualitative 
agreements between models are held better for smaller 
modulation depths of the potential Vq and for smaller 
pumping rates Pq. Eor a stronger pumping the oscilla- 
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FIG. 4: (Color online) Oscillation dynamics of the conden¬ 
sate calculated within the original model Q,([^. (a) The 

profile of the condensate density |T|^ of oscillating 7r-states 
for the potential contrasts Vq /h = 0.1 ps“^ and Pq = 18 
/im“^ps“^. (b) The phase profile of the dynamics shown in 
(a), (c) |T|^ profile of the meta-stable oscillations for the 

pump Pq = 30 /im“^ps“^. 
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FIG. 5: (Color online) Oscillation dynamics of the conden¬ 
sate calculated within the original model g,([§. (a) The 

|T|^ profile of the meta-stable oscillations for the potential 
contrasts Vo/h = 0.4 ps“^. (b) The phase profile of the dy¬ 
namics shown in (a), (c) and (d) The same as in (a) and 
(b), respectively, after a long time {t > 10000 ps). Other 
parameters: Pq = 18 /im“^ps“^. 


tions can become dynamically unstable and the system 
jumps spontaneously to the ground steady state with the 
minimal condensate energy [Fig. [^c)]. 

We also recalculated destabilization dynamics of the 
ground state for a stronger value of the potential con¬ 
trast Vo [Figs. (a),(b)]. The temporal beating between 
the ground eigenstate {jUgr) and the excited symmetric 
one {jJisym) can be clearly observed, which is in a full 















































































agreement with the results of the simplified calculations 
[see Eq. (16) and the Fig. (c) and (e)]. However, in 
accordance with the results of the numerical simulations 
of the original model 0, 0, this dynamical state is 
meta-stable. It means that after a long time (more than 
several tens of nanoseconds) the dynamics is transformed 
spontaneously [Figs.|^(c),(d)]. The analysis of the phase 
profile [Fig.[^(d)] shows that these hybrid oscillations ap¬ 
pear between Zero-state and 7r-states of the condensate. 
Apparently, the selective gain saturation associated with 
the reservoir dynamics is responsible for this particular 
choice of the condensate eigenstates involved into the os¬ 
cillations. We found that these types of oscillations are 
very robust and can develop even from the noisy initial 
conditions. 


VI. CONCLUSION 

In this work we considered nonlinear dynamics of co¬ 
herent exciton-polaritons within the periodical potential 
embedded into a planar microresonator driven by a ho¬ 
mogeneous incoherent pump. We restricted ourselves by 
weak-contrast lattices. Within this approach we devel¬ 
oped a simplified mean-field model for three spatial har¬ 
monics and found analytical expressions for the relevant 
eigenstates of the condensate. Detailed theoretical study 
of the condensate dynamics in the vicinity of the thresh¬ 
old Pth in the center and at the boundary of the Bril- 
louin zone was carried out. In particular, we demon¬ 
strated coherent persistent oscillations between different 
condensate eigenstates which are similar to those known 
for Josephson junctions. The numerical results have been 
supported by the analytical analysis. Also strong influ¬ 
ence of the dissipative effects and the incoherent reser¬ 


voir on the dynamics of coherent polaritons has been dis¬ 
cussed. 

Apparently, for strong nonlinear system confined in 
lattices, the influence of the exciton-polariton conden¬ 
sate nonlinearity on band structure becomes important 
and occurs with the increase of pumping rate Pq. As for 
atomic systems confined in a periodic potential it is well 
known that the inter-particle interactions can lead to the 
formation of loops (so-called “swallow tails”) both in the 
center and at the edge of the Brillouin zon^^H^. Such 
loops in the energy band structure have been interpreted 
as the indicators of superfluid condensate propertie^^ 
and hysteretic behavior of the condensate in respect to 
the variance of wave vectoi!^. More dramatic differences 
from the linear regime occur again in the center and at 
the edges of the Brillouin zone. Notably, such a system 
exhibits qualitatively new featrues of superfluid behavior 
of exciton-polaritons, affecting their energetic and dy¬ 
namical stabilit}®!. We hope to represent relevant results 
in the forthcoming paper. 
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